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ON THE PHASE-SPACE DISTRIBUTION OF BLOCH EIGENMODES FOR 

PERIODIC POINT SCATTERERS 


JORY GRIFFIN 

Abstract. Consider the 3-dimensional Laplacian with a potential described by point 
scatterers placed on the integer lattice. We prove that for Floquet-Bloch modes with 
fixed quasi-momenfum safisfying a cerfain Diophanfine condifion, fhere is a subse¬ 
quence of eigenvalues of posifive densify whose eigenfuncfions exhibit equidistri- 
bution in position space and localisation in momentum space. This result comple¬ 
ments the result of Ueberschar and Kurlberg lllTlI who show momenfum localisafion 
for zero quasi-momenfum in 2-dimensions and is fhe firsf resulf in fhis direcfion in 
3-dimensions. 


1. Introduction 

The phase space distribution of quantum eigenfunctions for large energies re¬ 
mains in general an unsolved problem - specifically, one would like to know whether 
the eigenfunctions of a given system exhibit equidistribution or some degree of lo¬ 
calisation (or indeed both). We are motivated by the physical problem concern¬ 
ing propagation through a cubic crystal lattice of scatterers. It is well known that 
when considering a scattering problem in which the wavelength is much larger 
than the radius of the scatterer, we can replace the scattering potential with a Dirac 
S point potential. This approach is perhaps most famously used in the Kronig- 
Peimey model [6]] which considers the one dimensional Schrodinger equation with 
a Dirac comb potential. A thorough treatment of models of this type, as well as 
higher dimensional analogs can be found in IH. Periodic problems of this sort can 
be tackled with Floquet-Bloch theory which allows us to reduce a periodic prob¬ 
lem in to a family of quasiperiodic problems on parametrised by their Bloch 
vector or quasimomentum k G T'^. 

For zero quasimomentum the problem of limiting phase space distributions has 
been studied in two dimensions by Rudnick and Ueberschar IfTTlI . and Ueberschar 
and Kurlberg HTTl |3, who showed that almost all eigenfunctions equidistribute in 
position space for all tori, and that in momentum space almost all eigenfunctions 
either equidistribute for square tori, or localise for tori with a diophantine ratio 
of side lengths. These results were partially generalised to three dimensions by 
Yesha |[T8llT9l who showed that for the cubic torus, all eigenfunctions equidistribute 
in position space, and that almost all eigenfunctions equidistribute in phase space. 
These results are further complimented by Kurlberg and Rosenzweig ||7|] who show 
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the existence of localisation in position representation in 2 dimensions and mo¬ 
mentum representation in both 2 and 3 dimensions. In this paper we generalise 
the results on the cubic torus to include nonzero quasimomentum, which destroys 
the high eigenvalue multiplicity and consequently, the equidistribution observed 
by Yesha. The proof follows a similar blueprint, but relies heavily on a result con¬ 
cerning the convergence of the two-point correlation function for inhomogeneous 
quadratic forms ||9l. 



(e) Aio4 « 100.11 (f) Aio5 w 100.13 

Figure 1. Six consecutive eigenfunction density plots on the plane 
(9,(p) showing the distribution of momenfum direcfions. We use fixed 
quasimomentum k = 


Problems of this type have been studied extensively in the Quantum Chaos liter¬ 
ature since Seba lil^ who considered a rectangular billiard with a point scatterer at 
some given point. The Seba billiard was constructed as an example of an infermedi- 
afe sysfem, meaning one fhaf is classically infegrable (the point scatterer affects only 
a zero measure set of trajectories) yet exhibits properties typical of chaotic systems 
||lHT2l,|T3Hlll. This is interesting in view of Shnirelman's theorem |0|l3[2O]] which 
states that for classically ergodic systems, a density one subsequence of eigenfunc- 
fions equidisfribufes in phase space, yef when the classical dynamics is integrable 
eigenfunctions tend to localise or scar. We are interested in the Schrodinger equa¬ 
tion on with potential described by point scatterers placed on 'Ird? which is 
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described by the formal operator 

(1) -A+c x: 

jelnZ? 

This operator is unitarily equivalent via a gauge transformation to a direct inte¬ 
gral over quasimomenta k. That is, we can instead consider a related quasiperiodic 
problem on the torus which is then realised via Von Neumaim self-adjoint exten¬ 
sion theory We first show that almost all of fhe eigenfunctions of fhis operafor 
equidisfribufe in position space. We then prove that there is a positive density se¬ 
quence of eigenfunctions which do nof equidisfribufe in momenfum space, specifi¬ 
cally we can find a subsequence fhaf partially localises in a given direction. 



(a) A]^4322 ~ 203.630 


(b) A23985 « 292.147 



Figure 2. A collection of non-consecufive eigenfunctions in momen¬ 
tum space with eigenvalue A showing partial localisation in the 
fixed direction (1, —1,0). We again use fixed quasimomenfum k = 

J_ X X) 

Vl' y/3' 
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2. Setup 


Consider the positive operator on = R^/ {2n'Z?) defined by 

( 2 ) + , 

The eigenfunctions of this operator are the complex exponentials 


(3) 


_ I _ 

( 271 ) 3/2 


with eigenvalue \^ + k\^, ^ G "Z?. Note here that when the components of k are 
irrational and linearly independent, then all such eigenvalues are distinct. This will 
be of major importance in the subsequent proofs. We will write J\f = {rij |; G N} to 
denote the ordered sequence of these eigenvalues. Equivalently, we could consider 
the standard Laplacian on T3 on functions that satisfy the quasiperiodic boundary 
conditions + 7 ) = for 7 in 2nZ?. In this case the eigenfunctions are 

proportional to the exponentials and again have corresponding eigenvalue 

1^ + /c|2 - it turns out that the first formulation is more convenient in our case. It 
is worth noting that this operator occurs naturally when considering the Laplacian 
on r 3 with some periodic potential. It is known that provided V{x + y) = V{x) 
for all 7 G 2tiI? then the operator on R3 given by —A + V{x) has a direct integral 
decomposition into operators on T3 of the form —A]^ + V{x). Full details of this 
procedure for a general operator can be found in IfTOl . We consider the perturbation 
of the operator — by a potential at a given point Xq G T3. We realise the 
perturbed operator 


(4) Hk — ~^k + 

via self-adjoint extension theory. Details of this calculation can be found in e.g. |l2l. 
The idea is that if we restrict our operator to functions vanishing at the point Xq, it 
should act like — Aj^. This operator is then positive symmetric but not self-adjoint, so 
we extend the domain of functions in such a way that self-adjointness is regained. 
If we define the restricted Laplacian, — Aq := — A with 

(5) Vq := C~(t3/{xo}), 


then the deficiency indices are (1,1) and the deficiency elements are the Green's 
functions, G±i(x, Xq), where we define G\ by 

(6) Ga(x,Xo) := (Afc +A) ^5{x-xq) = -^ ^ 

Throughout the paper we will also use gx = Ga/||Ga|| to denote the normalised 
Green's functions. There therefore exists a 1-parameter family of self-adjoint exten¬ 
sions parametrised by cp which we denote by The domains of these operators 
consist of functions / such that 

/(x) = C{cos{(p/2)^^^^^—^^+sm{(p/2)) +o(l) 


(7) 
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as X —> xq . The domain of ^ can be written 


(8) = |/i + cGi(-,xo) + ce^^G_i(-,xo) I h e Vq,c G C,^ G (-tt,7t)|, 

and the action of ^ is given by 

(9) -Afc^<^/= -Afc/i + ciGi(-,Xo)-ce^‘^iG_i(-,Xo). 

The new perturbed eigenvalues are given by solutions of the equation 


( 10 ) 

where 

( 11 ) 


E 

?ez3 


1 


1 ^ + ^1 


|^+?c| 2-A |^ + fc|4 + l 


C0= E 


= Cotan(^/2), 


e€Z3 


|^ + /C|4 + 1 


The set of perturbed eigenvalues will be denoted by A. 


3. Statement of Results 


We state the main results as two separate theorems, the first concerning pure 
position observables, the second concerning full phase space observables. To deal 
with phase space we first need to define quantisation. We follow the approach used 
in 11191 . Consider a classical symbol a G G“(S*T^), where S*T^ ~ x S^. We 
define the quantisation Op(ti) by 

(12) (Op(fl)/)(x) = e'(^'^^fl(x,^ + fc)/(0, 

?€Z3 


where we use the notation ^ 


X 

I?! 


and 1 ^ + is assumed to be nonzero for all 


We can then expand a in functions e^;^^(x,^) = Y) „j(^)e^(^'*), where is the 

(normalised) spherical harmonic of degree I and order m. Specifically we consider 
some finite polynomial P defined by 


(13) P{x,0= E E YL 

|f|<Ni 1<N2 \m\<l 


and claim that for all fl G G“(S*T^) there exist N\ and N 2 such that tor all (x,^) G 
and multi-indices oc with Ixl < 2 we have 


(14) 


- P{x,0)\ < £■ 


In light of this it suffices to prove our theorem only for these finite polynomials - the 
extension to a wider class of functions can be performed by expanding the function 
in a basis of these polynomials, truncating at some finite order, and controlling the 
error term (see I1T9I for details). We are now able to state the main results. Let A 
denote the sequence of perturbed eigenvalues. 


Definition 3.1. A vector k G is said to be Diophantine of type k if there exists a 
constant G such that for all m G G IN we have 


ifij C 

— > 
q q^ 


( 15 ) 


max 

i 
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The smallest possible value of k is k = 1 + ^. In this case k is called badly approx- 
imable. We now consider exclusively 3 dimensions, where k > 4/3. 

The first theorem concerns position space equidistribution and is proved in Sec¬ 
tion 5. 

Theorem 3.1. Fix cp ^ {—tz,tz). Assume the components of (1, k) are linearly independent 
over Q. Then, there is a density one subset A' C A such that for all observables a G C“(T^) 
loe have 

(16) lim {a{x)gx{x),gx{x)) = ^ ! a{x) dx 

A->oo OTT^ 

with A G A'. 

The second theorem concerns simultaneous equidistribution in position space and 
partial localisation in momentum space and is proved in Section 6 . 

Theorem 3.2. Fix (p G (—tt, tt). Let k be diophantine of type k G [4/3,2) and assume the 
components of {l,k) are linearly independent over Q. Then, for all e > 0 there is a subset 
Me C N of density at least 1 — e such that for all subsequences {M)neMe> there exists a 
further subsequence i^-nj)ieK sucfz that for all observables a G C°°(S*T^) we have 

(17) \im{Op{a{x,^))g.{x),gx {x)) = ^ / «(^,0 dxdy{^) 

where y has a positive proportion of its mass supported on a finite number of points. 

4. Truncation 

In order to consider only finite sums we define a fruncated Green's function. 
Define A(A, L) by 

(18) A(A,L) = {^ G : ||^ + fc|2-A| < L}, 
we then define the truncated Green's function by 

1 ^ Qi{^,x-Xo) 

(19) = 

and as before we denote by gA,L = Ga,l / || Ga,l || the corresponding normalised trun¬ 
cated Green's function. We want to show that for L = A"*^ for some d this truncation 
is a good approximation for large A. We first need a lower bound on the full Green's 
function. Define 

(20) J\f{x) = {n G A/" I n < x}. 

If the components of {l,k) are linearly independenf over Q then we know the as¬ 
ymptotic behaviour oi J\f (x) to be 

(21) N{x) = #J\f{x) = ^nx^^^ + 0(x®). 

It is conjectured that 6 = j + e for all e, and for k = 0, when counting with mul¬ 
tiplicities, the current best explicit bound due to Heath-Brown [|5l gives ^ = 52 + e 
for all e > 0. For our purposes if is required thaf 0 < 1, in fact we will show in the 
Appendix that we have 6 <i + £ independent of k. 
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Lemma 4.1. Let the components of {'l,k) be linearly independent over Q. Then, there is a 
density one subset of eigenvalues A' C A such that for A G A' 


( 22 ) 

Proof We have from ( [2l] > that 

1 


IGaII > 


(23) 




X 


N(x) 


X 


- 1/2 


Thus, since U]^ — > 0 we must have that for a subsequence of density one that 

(24) nfc+i - nfc < 

Let n]^_^_i > A > njt and we see 


(25) 


|GAf » E 


1 


> 


> 


1 


neAf 


(n - A)2 (nfc+i - A)2 (n^+i - n^,)^ 


> > A^ 


□ 


S'A — 

Ga 

Ga,l 

IIGaII 

Ga,l 


Ga 

Ga,l 


IIGaII 

IIGaII 


< 


< 2 


|Ga - Ga,lII 

IIGaII 

'|Ga — Ga,l 

IIGJI 


+ IIG 


Ga 

'A,l|| 


IG 


A,L| 


Lemma 4.2. Let L = X ^, then ||yA,L “ S'aII —> 0 as A —^ oo with A G A’. 
Proof First we see that 

(26) ||g, 

(27) 

(28) 

(29) 

Then we have 

(30) 


Ga,l Ga,l 
1 


1 


IG 


A I 


G 


A,L| 


|GA-GA,L|r« E 


|f+fc|2-A|>L 


|^ + ;c| 2 -a) 2 - 


We evaluate the lattice sum via Abel summation, which tells us that for a smooth 
function f we have 


(31) E f(\i + T) = N{nB)f{nB)-N(nA)f(nA+i) 

WA<l?+fc3<WB 


rtlB 


'nA+1 


fft)N{t) dt. 


mB 


Integrating by parts we see 

(32) E /(l? + ':|")=2ir/'" 

«A<|?+fc3<WB 

(33) +0{nlf{nB) -n^Afi^A+i)) + 0{ f \f'{t)\t^ dt). 

JnA+i 
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Applying this to f{n) = with = hq and ng < — L < ns+i we see 


_ - _ = 271 

(n - A)2 ^ (n - A)2 

We can bound the integral by 

rns n'^/2 


-Wb ^ 1/2 


dn + O 


(ng - A)^ 


r-WB 


(A — ny 


dn . 


- -^ dn < A^^^ 

(n — A)^ 


nWB 1 


{n-xy 


<A!<Y 

- L - 


Similarly we see 


and also 


_ _ < _ 

(nB-A)2 - L2' 


/„! (A-n)3 L2- 


Now repeating this procedure with < X + L < and np = oo we obtain 
(39) 


_ - _= 2n 

n^-L (” “ '^^A+i (n - xy 

For the first integral we write 


noo ^^/2 


dn + 0 


{ha— 


(A — n)3 


dn . 


rco ^3/2 


(n — A)^ 


dn = 


/■“ (s +A)^ 

/w^+i-A s2 


^ f^{s + xy/ 
- Jl s2 

(40) « — « 

For the second term we have immediately 


(s + A)^ 


/A s" 


- — - — . 

(n^-A )2 L 2 


For the third term we see 


r-OO 


(A — n)3 


dn = 


(s + A)^ 


/■“ s + A 


'wa+ 1 -A S 


i^ds 

A 


Putting all of this together we see 


Ga-Ga,l|| ^72' 
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and hence, using Lemma 4.1 that for the normalised Green's functions 


(44) 


\gA,L-gA\\ < 


L 


1 

which tends to 0 for all S — e. 


Corollary 4.1. Define as above with L = A ^ and 0 < — e then 

(45) I {OY>{e^,i,m)gx,Lrgx,L) - {Ov{^Q,m)g\rg\) I ^ 0 . 


□ 


Proof. We have 

(46) I (Op(e^,/^^)gA,L/gA,L) - (Op(e^;,m)gA/gA) I 

(47) < I (Op(e^;,^)gA,L,gA,L - ^ a ) I + I (Op(e^,,,^)(gA - gA,L),gA)\- 

Taking each ferm and using Cauchy-Schwarz gives 

(48) |(Op(e^,z,^)gA,L/gA,L) - (Op(e^,/^;„)gA/gA)| < II Op(e^;,^) || ||gA - ^aaII ^ 0. 

□ 


5. Equidistribution in Position Space 
The following proposition is key fo the proof. 

Proposition 5.1. Fix ^ 7 ^ 0, Z G N and \m\ < 1. Let L = for some <5 > 0. Let the 
components of (1, k) be linearly independent over Q. Then, for A sufficiently large we have 

(49) {0^{^^,l,m)gA,L>g\,L) = 


In order to prove this we first need a lemma. 

Lemma 5.1. Let the components of {k,l) be linearly independent over Q and fix ^ G 
nonzero. Then, there exists some e > 0 such that for all ^ G Z^ we have 

(50) \2{^ + k,0 + \ff\>e. 


Proof. We have that 

m + U) + \l\^\ > ||2{*:,f>ll 

where || • || represents the distance to the nearest integer. Since we assumed the 
components of {k,l) were linearly independent, this is bounded away from zero. 
Now choose e = II (fc,^) II. □ 
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Proof of Proposition \5. 1 \ First write 
I (Op(e^,/,m)GA,L/ Ga,l) I 


(51) 


(52) 


1 


64n^\\G\ IP 


E 


?€A(A,L) 


ei(?A-xo) 

I? + )cF-a 




!?€A(A,L) 


„i()7,x-xo) 

^ \ 

+ fcp — A' 


647r^||GA IP 


/T3 


E 




^,n^A{KL) 


(l^ + ^p-^)(l7 + ^p- 




Integrating over x leaves only the terms where r] = f, P f. However, note that by 
Lemma 5.1 for ^ G A(A, L), 

(53) ||^ + ^ + ;c| 2-A| = ||^ + fc|2_A + 2(^ + ;c,0 + |^|"| »e 

so ^ ^ ^ A{\, L) for A sufficiently large. Thus the integral in ( [5^ vanishes. □ 

We are now able to show equidistribution for position space observables. 

Proof of Theorem 3.1 Let A G A', and let a G C°°(T^). The operator Op(fl) is then just 
given by multiplication by a. We consider a to be some finite polynomial 


(54) 


i(^)= E 

T\<N 


and see from Proposition |5.1| that 

(55) {a{^)gx,Lrg\,L) {K^)gKLrg\,L 

The result then follows from Corollary |4.1[ 


/t3 


\gA,L{^) 


2 dx 


8n^ 


□ 


6. Localisation in Momentum Space 

Throughout this section we will assume k is diophantine of type k < 2. Let a be 
defined by 

(56) a{x,^)= d{f,l,m)e^^i^^{x,f,) 

\l,\<Ni,l<N2,\m\<l 


where d{f,l,m) is given by 

(57) d{f,lm) = ^ J^^a{x,^)e-^<^'^'>Yl^{l)dxda{f,). 
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We thus have that 


(58) (Op(fl)gA,L,gA,L) ^ (ER(0,/,m)Op(eo,/,m)gA,L,gA,L) 

l,m 


= iiGAir"7i^E E »(0T, 


6471^ 


m] 


l,m ^eA{\,L) 


YiA^ + k) 
A + A-Ay 


|Ga|| ^ E L a{xAdx 

e6A^,L)-^S*T3 (|^+fc|2_A)^ 


Thus the component of the spectral measure for each fixed |^ + fcp = m on x 
consisfs of Leb x 5-^^. The full (uimormalised) spectral measure is thus a weighted 

sum of a growing number of 5 masses fhaf become dense on S^. We aim fo show 
that for a positive density subsequence of A, the tails of this sum can be bounded 
uniformly in A such that a positive proportion of ifs densify will be supporfed on a 
finite number of poinfs. A key feafure of the proof is the work by Marklof 0. We 
now recapitulate the necessary results in dimension 3. 


Definition 6.1. Let M = {uj | ; G N} be defined as before. Let 1 / 71 , 1/72 G »S(R+) be 
Schwartz functions, and let h G Co(IR) be compactly supported, and h its Fourier 
transform. Define fhe generalised pair correlation function by 

(59) R{,pi,,p2,h,r) = E ^>1 {”f) li(VT(n,.-«;)). 

%i 


Theorem 6.1 (See Theorem 2.5 in [|9l for the statement in full generalify). Let k be 
dioTphantine of type k < 2, and that assume the components of (k, 1) are linearly independent 
over Q. Then 

,i/72,/i,T) = J 



poo 

h(s) ds / \pi{r) xp 2 {r) r dr. 
Jo 


We now proceed with the proof. 

Lemma 6.1. For G > l,we have that #{n; G J\f{T) : n^+i — n,- > G/ y/fffG} < T^^^/G. 


Proof We see thaf 
(60) 


E Ai+iA+i - ni) < E All - 

ni<T ni<T 


(61) < 

Thus by Chebyshev's inequality we see 

(62) #{ni < T : Si = — ni>G/y/nffi} < G. 


Lemma 6.2. Given D > 0, E > 1, 

(63) #{n G AA(r) : \M{T) n [n - n + 


> E + 1} < 


Dr3/2 


□ 


E 
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Proof. We have that 

(64) (|AA(r)n[n-^,n + ^]|-l) 

neAfiT) 

(65) = #{n,m ^ J\f{T) : m n, ^yn\n — m\ < D} 

( 66 ) #{n,m G J\f(T)\J\f{T/2) : m n, Vf\n — m\ < D}. 

Since we assumed k was diophantine, we may apply Theorem [O with ipi = ip 2 the 
indicator function of [1/2,1], and h the indicator function of [—D, Dj. This gives us 
the asymptotics 

(67) #{n, m G Af {T)\Af (T / 2) : m ^ n, \/f\n - m| < D} - Stt^DT^/^. 

Again by Chebyshev's inequality we conclude 


( 68 ) 


#{n e V(T) : |V(r) n (« - i n + i] > E + 1 } 


D 




DT^/2 


□ 


Lemma 6.3. For all A > 1 
(69) 


E 

n,mEAf{T) 

m\n—m\>A 


m(n — m] 




p3/2 

7^’ 


Proof. We first define 

(70) M{k) :=\{n eM e[Kk + l]}\. 

Then we deduce an bound on M{k) by 

(71) M{kf= Yu ^[k,k + \\}\ 

fc<T3/2 fc<T3/2 

(72) < I {m, n G A/": + 1, e [-1,1]} | 


which again by Theorem 6.1 gives us 
(73) Y ^ 

fc<T3/2 

Nofe fhat we can write 

_ 

(y/m\n — + — i^l) > 


m n — m =—= 


7m + 7n 


m + 7n 


|n3/2_^3/2| 


and also that 7m\n — m\ < Hence, we can bound the sum in ( [69l ) by 

1 - (i + Gi: 


E 




E 


n,meAf{T) n,meM{T) 

y/m\n—m\>A \n^^^—nP^^\>A 


W\1 

_ m ‘ 

n^/2 _ j^3/2^2 

(1 + 


[T3/2J 

E E 

k=A n,m€Af{T) 


n-\2 
m ' 


3/2 _ 


j^3/2^2 


( 74 ) 
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Now, when m > n we can immediately conclude 

[T3/2J 

E 


E 


k=A 


(75) 


n,mEAf{T) 


{n 


(I _L /^')2 -I 


3/2 _ 


m 


3/2y. 


k=A 


n^^^-m^/^\e[k,k+l] 


k — jA 

By Cauchy-Schwarz we may bound this above by 

1 7 - 3/2 

E 


<4 E ^ E M{m){M{m + k) + M{m + k + l)). 


k=A 


When m < n, we see that 

(76) -m^'^ ^[k,k + l] 


A 


n\3/2 fc + 1 
-I <1 + 


m) 


m 


3/2 ’ 


We know that m is bounded away from zero, say m > C, then we must have 
(77) 


mJ 

Repeating the previous argument in this regime yields 


[T3/2J 

E E 

k=A n,mEJ'/{T) 


(1 + 


W\2 


[r3/2j 


1 


(j^3/2 _ ffi3/2y 


< E y^\{m,neAf{T):{n^/^-m^/^)e[k,k + l]} 


k=A 


n^^^-m^/^\e[k,k+l] 


(78) 


(79) 


LT3/2J 

< E 


4/3 E M{m){M{m + k) + M{m + k + l)) 
k=A ^ m<T^/^ 


1 7-3/2 

7-3/2 y- _L ^ 1 _ 

^ j^A/3 ^ ^ 1 / 3 ' 


fc=A 


□ 


where the final line follows from Cauchy-Schwarz as before. 

We are now ready to prove the second main theorem. 

Proof of Theorem |3.2[ Define J\f' as follows, first remove from J\f all points m whose 
nearest left neighbour is further than G/ \/m, by Lemma | 6 .l| we are left with a 
subsequence of density at least 1 — 1/G. Now choose D and fix E large enough 
such that 

7-3/2 

(80) '' ‘ ^ '■ ^ ^ - 


\{meM{T) : \M{T) n [^- > £ + 1} 


< 


G 


which is possible by Lemma 6.2 Removing these points leaves us with a subse¬ 
quence of density at least 1 — 2/G. Finally, by Lemma 6.3 and Chebyshev's in¬ 
equality we choose F large enough such that 

1 t 3/2 

(81) \{meAf{T): ^ 

neM{T) ' ' 


G 
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Removing these points leaves us with a subsequence of density at least 1 — 3/G. 
Thus if we consider pure momentum observables and for m G J\f' denote by }irn the 
delta measure on the point corresponding to the direction ^ + k with |^ + fcp = m, 
we see that the unnormalised measure associated to G a is 

/Vm 

_ V-m _I y V-n y j^n 

{m-XmY {n - XniV {n - X^nV' 

n I I D I I D 

0<\n—m\<—r= \n—m\>—r= 

' ym ' ' ym 

We know that the mass of the first term is S> m/ G^, the mass of the second sum 
has at most E terms, and the mass of the third is bounded above by F m. Thus the 
normalised measure will have a positive proportion of its mass on a finite number of 
points. The theorem then follows from compactness of S*T^ and by setting e = 3/G 
and defining Me by j G Me ny G J\f' where ny is the ordered unperturbed 

eigenvalue. □ 


Appendix A. 

Proposition A.l. Let S{R) = #{|^ + fc| < R | ^ G denote the number of shifted 
lattice points inside a ball of radius R. Then we have that 

(83) S(R) = ^nR^ + 0{R^^^+^). 

Proof We bound the quantity S{R) above and below by sums over the indicator 
function of a shifted ball convolved with some smooth bump function with smooth¬ 
ing parameter 5. We can then employ Poisson summation and tune 5 in such a way 
that the error terms vanish. Let denote the ball of radius R centred at k, and 

write tps{x) = S~^xp{x/S) where if) is some smooth function with compact support 
in Bo(l) normalised such that ip{0) = 1. Define Sj(R) to be the smoothed sum 

(84) Ss{R) = XBXR)*fs{x). 

Note that we have S{R — S) < Ss{R) < S{R + S). By Poisson summation we see 

(85) E ?Cbxr)*M^) = E ^bxr)(^)M^)- 

xEl? 

Computing the term ^ = 0 yields 

r 4 o 

(86) J^^Xbxr){x) dx = -ttRT 
For ^ 7 ^ 0, the Fourier coefficients Xb^r) (^) given by 

(87) XbXR)( 0 = 27T^j^j3 (sin(27rR|^|) -27rR|^|cos(27rR|^|)). 

We also have that 

[ 5-^ip{x/5)e-^^^^^'^Ux = [ 5-^ip{x/S){An^\^f)-\-A)e-^^^^^'^Ux 

( 88 ) = [ r^e-^^^^’''^^(-A)ip(x/d)dx 
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We thus have that \^s{0 I mm{l, (^1^1) < (^1^1) Plugging these 

asymptotics in to the sum gives 

(89) E Xb,(r) * E 

x6Z3\{0} ?€Z3\{0} 


We thus have that 

(90) S{R) < Ss{R + S) = ^n{R + + 0{Rr^^+^^) 

= ^nR^ + 0{R^S + 

and similarly that 

(91) S{R) > Ss{R -S) = ^nR^ + 0{R^5 + R^-(i+^)). 
Setting S = R~^^^ yields the result. 


□ 
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